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Abstract: Mathematical modeling of many physical systems leads to nonlinear evolution equations because most 
physical systems are inherently nonlinear in nature. The investigation of traveling wave solutions of nonlinear partial 
differential equations (NPDEs) plays a significant role in the study of nonlinear physical phenomena. In this article, 
we construct the traveling wave solutions of modified KDV-ZK equation and viscous Burgers equation by using an 
enhanced (6 7G) -expansion method. A number of traveling wave solutions in terms of unknown parameters are 
obtained. Derived traveling wave solutions exhibit solitary waves when special values are given to its unknown 
parameters. 
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Background 

Engineers, physicists, and mathematicians have always 
shown their incessant interest in studying nonlinear 
problems related to numerous scientific applications, 
such as fluid dynamics, high-energy physics, plasma 
physics, elastic media, optical fibers, biomathematics, 
chemical kinematics, chemical physics and geochemistry. 
Many young scientists have also shown their increased 
interest in last two decades because of plausible develop- 
ment of nonlinear science during this period of time. In 
order to understand the behavior of a nonlinear phe- 
nomenon we need to solve the nonlinear equation/set of 
equations describing that phenomenon, is often very much 
challenging. There are so many approaches developed 
over years to analyze/solve such system of nonlinear 
equations, most of them are based on some assumptions, 
and hence approximations. Though perturbation methods, 
like other nonlinear analysis techniques, have their own 
limitations, are most useful methods among all these 
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approaches so far. Using perturbation method, to achieve 
the ideal results an appropriate choice of small parameter 
has to be made efficiently, otherwise, a fatal error in re- 
sults may occur. The perturbation methods are not even 
applicable to many nonlinear equations because of not 
having small parameter, which is the principal assumption 
that has to be met by an equation in order to apply per- 
turbation method. Furthermore, the approximate solu- 
tions obtained using perturbation methods are valid only 
for the small values of the parameters (Ghorbani and 
Saberi-Nadjafi 2007; Mohiud-Din 2007; Mohyud-Din and 
Noor, 2009). The investigation of exact traveling wave so- 
lutions to these nonlinear equations (NPDEs) have also 
been observed as a field of great interest to many mathe- 
maticians and physicists because of its significant role in 
understanding the behavior of nonlinear physical phenom- 
ena. As a result, numerous techniques of obtaining travel- 
ing wave solutions have been developed over last three 
decades, such as, the Hirota's bilinear transformation 
method (Hirota 1973, 1981), the modified simple equation 
method (Jawad et al. 2010; Khan and Akbar 2013a; Ahmed 
et al. 2013; Zayed and Hoda 2013; Zayed and Arnous 
2012), the tanh-function method (Wazwaz 2005; Parkes 
and Duffy 1996), the Exp-function method (He and Wu 
2006; Akbar and Ali 2011b; Bekir and Boz 2008; Xu et al. 
2009), the Jacobi elliptic function method (Ali 2011), the 
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(G VG) -expansion method (Akbar et al. 2012a, 2012b; 
Akbar and Ali 2011a; Wang et al. 2008; Shehata 2010; Koll 
and Tabi 2011; Naher et al. 2013; Zayed 2009, 2010; Asian 
2010; Bekir and Aksoy 2012), the homotopy perturbation 
method (Mohiud-Din 2007; Mohyud-Din and Noor 
2009), transformed rational function method (Ma and Jyh 
2009; Ma and Fuchssteiner 1996), multiple exp-function 
method (Ma et al. 2010; Ma and Zhu 2012), gene- 
ralize Hirota bilinear method (Ma 2013), enhanced 
(G VG) -expansion method (Khan and Akbar 2013b), The 
Sine-Cosine method (Bibi and Mohyud-Din 2013), the 
first integral method (Tascan and Bekir 2010; Feng 2002), 
the ansatz method (Hu 2001a, 2001b) and many others. 

The present article is devoted to construct the exact 
solutions for modified KDV-ZK equation and viscous 
Burgers equation using a relatively new technique, 
named, enhanced (G VG) -Expansion method. The rest 
of the paper is organized as follows. Details explanation 
of enhanced (G VG) -expansion method has been pre- 
sented in the next section. The obtained solutions of 
modified KDV-ZK equation and viscous Burgers equa- 
tion using this method are presented in Section 3. In 
Section 4, we have presented some graphs of obtained 
family of solutions for some particular values of the un- 
known parameters and the final conclusions are shown 
in Section 5. 

Enhanced (G VG) -expansion method 

In this section we describe in details the enhanced 
(G VG) -expansion method for finding traveling wave 
solutions of nonlinear evolution equations. Any nonlin- 
ear evolution equation in two independent variables x 
and t can be expressed in following form: 



R(u, U t , U X7 U t t 1 Uxx, u xt , ) — 0, 



(2.1) 



where u(£) = u{x, t) is an unknown function, R is a poly- 
nomial of u{x, t) and its partial derivatives in which the 
highest order derivatives and nonlinear terms are in- 
volved. The following steps are involved in finding the 
solution of nonlinear Equation (2.1) using this method 
(Khan and Akbar 2013b): 

Step 1: The given PDE (2.1) can be transformed into 
ODE using the transformation <f = x ± Wt, where W is 
the speed of traveling wave such that Ws R - {0}. 

The traveling wave transformation permits us to re- 
duce Eq. (2.1) to the following ODE: 



Step 2: We now suppose that Eq. (2.2) has a general 
solution of the form 



"(0 = £ I 



MG'/G)' 
(1 + X(G'/G)Y' 



h{G'/Gf 



(G'/cy 



(2.3) 



subject to the condition that G = G{£) satisfy the equation 

G" + AG = 0 (2.4) 

where a it b i{-n < i < n; n e N) and A are constants to be 
determined, provided that a = ± 1 and \i * 0. 

Step 3: The positive integer n can be determined by 
balancing the highest order derivatives to the highest 
order nonlinear terms appear in Eq. (2.1) or in Eq. (2.2). 
More precisely, we define the degree of w(£) as D(u{£)) = 
n which gives rise to the degree of other expression as 
follows: 



D 



d? 



q 7 D[ uP 



np + s(n + q) 

(2.5) 



Step 4: We substitute Eq. (2.3) into Eq. (2.2) and use 
Eq. (2.4). We then collect all the coefficients of (G VGj 



and 



{G'/Gl^ 



(G'/G)' 



together. Since Eq. (2.3) is 



a solution of Eq. (2.2), we can set each of the coefficient 
equal to zero which leads to a system of algebraic equa- 
tions in terms of a L , bi{-n <i<n;n& N), A, and W. One 
can solves easily these system equations using Maple. 
Step 5: For fi < 0 general solution of Eq. (2.4) gives 



G 
G 



and 



G 
G 



tanh(j4 + s/^C) 



-^coth(^4 + yf^C) 



And for pi > 0, we get 

G 



and 



G 



^tan^-y^) 



^cot^ + y/jit) 



(2.6) 



(2.7) 



(2.8) 



(2.9) 



R(u,u',u", ) = 0, 



(2.2) 



where R is a polynomial in u(Q and its derivatives, where 
M'(f) =f,w"(0 =|f, and so on. 



where A is an arbitrary constant. Finally, we can con- 
struct a number of families of traveling wave solutions of 
Eq. (2.1) by substituting the values of a- u bi{-n < i<n;ns 
N), A, and W (obtained in Step 3) and using Eqs. (2.6) - 
(2.9) into Eq. (2.3). 
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Applications of enhanced (G VG) -expansion method to 
modified KDV-ZK equation and viscous Burgers equation 
Modified KDV-ZK equation 

In this current sub-section, we apply enhanced (G VG) - 
expansion method to solve the modified KDV-ZK equa- 
tion of the form, 



Set 4: 2? = 0, W = 6fi,A = a. x = +-_ 



: 2d a , — + Jb_ 

do = ao, a\ = 0, b.\ =0, bo = 0, b\ = 0. 
Set 5: R = 0, W = -3p, X = X, = 0, a 0 = 0, 

a x = 0, b.i = 0, b 0 = 0, h = ± 73 
Set 6: A' 0. IV 



-3^, A = A, a_i = 0, «o = 0, 



a\ = 0, b.\ = 0, 



I6fi 
Plod ' 



bi = 0. 



u t + du u x + u x 



0 



(3.1) 



where d is a nonzero constant. 

The traveling wave transformation equation = u 
{x,t), £=x+y+z- Wt transform Eq. (3.1) into the fol- 
lowing ordinary differential equation: 



-Wu' + du 2 u' 



0 



Integrating Eq. (3.2) with respect to we obtain 
3u"-Wu + ldu 3 + R = 0 



(3.2) 



(3.3) 



where R is a constant of integration. Following the 
process as described in step 3 (Section 2), balance be- 
tween the highest-ordered derivative term u" and non- 
linear term u 3 of Eq. (3.3) provides n = 1. 
For n = 1 Eq. (2.3) takes the following form: 



u(Q = a 0 + 



a x {G'/G) 



1 + A(G'/G) 
MG'/G)" 1 



a. x (l + X{G'/G)) 
(G'/G) 



a\ 1 



(G'/G) 2 



b x 



\ 



a\ 1 



(G'/G) 2 



MG'/G)" 



CT I 1 



(G'/G) 2 



(3.4) 



where G = G(0 satisfies Eq. (2.4). Substituting Eq. (3.4) 
into Eq. (3.3) and using Eq. (2.4), we get a polynomial 



1 



G'/GY 



in (G'/GY and {G'/G)'^a^ 
coefficient of (G'/G/ and {Gt/GfJo^ 



Setting the 



(G'/Gy 



equal 



to zero, we obtain a system containing a large number 
of algebraic equations in terms of unknown coefficients. 
We have solved this system of equations using Maple 
13 and obtained the following set of solutions: 

Set 1:2? = 0, W = -12/a, X = 0, a. x = a 0 = 0, 



/2d ' 



= 0, 6o = 0, bi = 0. 



Set 2: R = 0, W = 6fi, A = 0, a„i = 0, a 0 = 0, 



/2d' 



Set 3: i? = 0, W : 



b. x = 0, & 0 = 0, bi = 0. 

A = 0,a^i = 0, ao = 0, 



. 2c/ 



fc.i = 0, 6 0 = 0, fcj 



2 - 



Corresponding to each set of solution we get two dif- 
ferent families of traveling wave solutions of Eq. (3.1) ac- 
cording as fi < 0 and fi > 0. In doing so, we have obtained 
twelve families of traveling wave solutions while each 
family comprises of two solutions. 

First we represent the families of hyperbolic solutions 
corresponding to fi < 0: 

Family 1 : = ±3 ^ ( tanh (A + s /=j%) - i coth (a + y^if) ^ , 

« 2 (f) = ±3 \J^( coth(A + </=!%)— tanh(A + y^f))' 

where f=x + y + z+ 12^f. 



Family 2: « 3 (£) = ±3^/^ tanh (4 



coth (/I 



«4 (0 = ±3^ 
where <f=x + y + z- 6/^£. 

Family 3: « 5 (f) = ±3W^( tanh(A + y^if) + Isech(A + yS<0), 
u 6 (f) = +3y/^(f coth(A + y/^/Zf) + cosec/)(A + y^f)), 

where £ = ;t + y + z- § ^t. 

Family 4: M 7 (f) = 2a 0 +3y , '-^coth(A + , 

w 8 (0 = 2a 0 +3^tanh(A + V^), 

where ^ = ^ + ^ + 2- 6fit. 

[2~ _ _ 

Family 5: "9(f) = +I3J — cosech(A + y^f) coth(A + y^f), 

/~2~ _ _ 

"10(f) = +3J—sech(A + y/=jK) tanh(A + y^f), 

where £=x + y + z + 3fit. 



Family 6: Un{£) = +3\ -^-cosech(A 
yd 

u n (n = +31 J^sechtA- 
V d 

where £,=x + y + z + 3^t 

And, the families of plane periodic solutions 

corresponding to [i > 0 are: 

Family 7: «i 3 (0 = +/3y^^tan(A-^) +i cot(A-0Sf)J , 

« 14 (0 = +^J^( + v^«f) +^ tan ( A + vK)J . 
where f=x + y + z+ 12^t 

Family 8: Wi 5 (<9 = T/Sy^tan^-y^f), 
+/3J^cot(A+ vK), 



"16(f) 

where f=x + y + z- 
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Family 9: u 17 (Q = +/3^(tan(A-^) + sec^-^)), 



= +3 y^( cot (^ + v^) + cosec(^l + rftf)), 



where £ = x + j + z- 1 



Family 10: = 2a 0 +I3J-jCot(A- y /fi^), 

«2o(0 — 2a 0 +I3^j^tan(A + y^ - ), 
where £=x+j+z- 6^t. 

/~2~ 

Family 11: "21(f) = ±13 J— cosec(A-^Sf) cot(A-0«£), 
K 22 (0 = ±I3^sec(A + ^ ) tan(A + , 



where £ = x + j + z + 3^£. 



Family 12: W23(£) = +/3W — cosec^-y^f), 
i*m(0 =T/3W^sec(A + ^), 



Remark 1: All the obtained solutions have been 
checked with Maple by putting them back into the ori- 
ginal equations and found correct. 

lA'scous Burgers equation 

In this sub-section, we will apply enhanced (G '/G) - 
expansion method to solve the viscous Burgers equation 
of the form, 



u t ± uu x -vu x 



0 



(3.5) 



where v is the viscosity coefficient. Burgers equation is a 
model for nonlinear wave propagation, especially in fluid 
mechanics. It occurs in various areas of applied mathem- 
atics, such as modeling of gas dynamics and traffic flow. 

The traveling wave transformation equation = u 
(x, t), f = x - Wt transform Eq. (3.5) into the following 
ordinary differential equation: 



where £ = x + y + z + 3fit. 



-Wu' 



(3.6) 
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Integrating Eq. (3.6) with respect to f setting constant 
of integration to zero, we obtain 



2vu' + 2Wu-u 2 = 0 



(3.7) 



Following the process as described in step 3 (Section 2), 
balance between the highest-ordered derivative term u 
and nonlinear term u of Eq. (3.7) provides n = l. 
For n = 1 Eq. (2.3) takes the following form: 



a x (G'/G) 0_i (1 + A(G'/G)) 



1+\(G'/G) 
+b 0 (G'/G)- 1 



(G'/G) 



a\ 1 



(G'/G) 2 



a\ 1 



(G'/G) 2 



(3.8) 



+b- 1 (G , /Gy 



\ 



o- 1 



(G'/G) 2 



where G=G(f) satisfies Eq.(2.4). Substituting Eq. (3.8) 
into Eq. (3.7) and using Eq. (2.4), we get a polynomial 



coefficient of (G'/G) 7 ' and (G^G)' J + 

equal to zero, we obtain a system containing a large 
number of algebraic equations in terms of unknown coef- 
ficients. We have solved this system of equations using 
Maple 13 and obtained the following set of solutions: 

Set 1: W = ± Vy^JI, A = 0, «_i = 0, «o = ivy^, 
a\ = -v, b-i = 0, bo = 0, b\ = ±V\ 



Set 2: = +2v^/— jS, A = , a-\ = 0, a§ = ag. 

2vpi 



2vpi 



0. 



0, 



0. 



Set 3: W = ±2v. v / z ^, A = 0, a_i = 2v^, 0 O = ±2v 
y 3 /^, a.\ = 0, &_i = 0, bo = 0, &i = 0. 
Set A: W = ±4v v / = /^,A = 0,a_i = 2v,«,«o = i'lv^/ 3 ^, 
a x = 0,^_i = 0,^o = 0,*i = 0. 
Set 5: W = ±2v^, A = - ( ^Zf^ ) > «-i = 



^0? ^1 



2v^ 

0,fc 0 > 0,6; 

0 



0, fc-i 

Set 6: W=±vy=^, A = -(^ 
<2o = 0o, 0i = 0, 6-i = 0, bo = 



-- 0. 

0-1 = v/*, 
6i = 0. 



in (G 7G) j and {G'/G)' Ja(l + . Setting the where £ = x+v^ 



According to the parallel course of action the families 
of hyperbolic solutions corresponding to fi < 0 are: 

Familv 1 • = v^/=ji(l- UnhU + y/=j% ) +1 sech (A + V^f)), 

' u 2 (Q = v^j2(l-coth(A + ^)Tcsch(A + y/=t£)), 



■fit. 
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Figure 4 Periodic wave profile of mKDV-ZK equation. 
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Family 2: U3 (Q = 2a 0 vfi 

"4(f) = IUqV^I 



1 ± tanh(A + J=tf) 
2vfi + (fl 0> /-^+2vf<) tanh(A + ^jiCj 

1± coth(.4 + y/=f%) 
1v\i + («o-\/-^+2v^) coth(A + 



= x ± 2Vyf zr jlt 
M 5 (£) = ±2v 



^(l+coth(^4 
^(l+tanh(y4 



« 6 (0 = ±2v 

= X+lVyJ^flt. 

Uj{£,) = ±2v 
w 8 (£) = ±2v % 

Ug{^) = ±2v\ 

«io(0 = ±2v 

«n(0 = v^/^fl-cothfA + y^^TcscA(A+ x/^f) 
"12(f) = v,/=^(l-tanh(^4 + V / Z /Jf) Wsecfe(A + ^/t^ 



^(2+cothiA -| 
^(2+tanh(A 4. 

^(l+coth(/t 4 
^(lTtanh^ 



where £ = 
Family 3: 

where £ = 
Family 4: 

where £ = 
Family 5: 

where £ = 
Family 6: 

where £ = 



Similarly, we can write down the other families of 
exact trigonometric function solutions corresponding to 
fi > 0 which are omitted for convenience. 

Remark 2: All the obtained solutions have been 
checked with Maple by putting them back into the ori- 
ginal equations and found correct. 

Graphical representation of solitary waves 
(obtained from above family of traveling waves) 
Modified KDV-ZK Equation 

We have obtained total twenty four solution profiles of 
traveling waves in terms of some unknown parameters; 
subdivided (above) into twelve families according as 
negative or positive values of p. All of these solutions 
are combination of hyperbolic functions or trigonomet- 
ric functions but Family-4 and Family-10 which are 
combination of algebraic functions, and hyperbolic func- 
tions or trigonometric functions. Solitary waves can be 
obtained from each traveling wave solution by setting 
particular values to its unknown parameters. If we put 
«o = 0, then Family 4 and Family 10 coincide with Fam- 
ily 2 and Family 8 respectively. In this section, we have 
presented some graphs of solitary waves constructed by 
taking suitable values of involved unknown parameters. 

• For the values of ^ = - 1, d = - 2, A = y = z = 0 
within the interval - 5 < x, t < 5 Figure 1 is kink 
wave (only shows the shape of Us(Q of mKDV-ZK 
equation). 

• For the values of ^ = - 1, d = I, A = y = z = 0 
within the interval - 5 < x, t < 5 Figure 2 is soliton 
wave (only shows the shape of of mKDV-ZK 
equation). 

• For the values of fi = 0.50, d = - I, A = 3, y = z = 0 
within the interval - 5 < x, t < 5, Figure 3 is a 



periodic wave (only shows the shape of Wis(£) of 
mKDV-ZK equation). 
• For the values of = 1, d = - 2, A = y = z = 0 
within the interval - 5 < x, t < 5, m 2 4(£) (solution of 
mKDV-ZK equation) shows the shape of periodic 
wave represented in Figure 4. 

Viscous Burgers equation 

Now we will discuss some of the obtained results of vis- 
cous Burgers equation and their graphical representations. 

For the values of A = 0, (i = - 2, a = 1 and v = - 2 
within the interval - 3 < x, t < 3 the solution of viscous 
Burgers equation Mi(£) shows Kink wave which is repre- 
sented in Figure 5 and For the values of «o = 0, A = 0, 
ft = - 1, a = - 1 and v = - 0.50 within the interval - 3 < x, 
t < 3 the solution of viscous Burgers equation u^ity shows 
singular kink wave which is represented in Figure 6. Some 
graphical representations of viscous Burgers equation are 
given below: 

Conclusions 

In this article, we obtain a large number of exact travel- 
ing wave solutions including solitary wave solutions for 
modified KDV-ZK equation and viscous Burgers equa- 
tion through enhanced (G 7G))-expansion method. 
Using this method, we have found some extra family of 
solutions, is the best thing about this method. The ob- 
tained solutions suggest that enhanced (G 7G))-expan- 
sion method can be used as useful mathematical tool for 
solving nonlinear evolution equations arises in the arena 
of mathematical physics, engineering sciences and ap- 
plied mathematics. 
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